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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this 
to  metric  (SI)  units  as  follows: 

_ Multiply _  _ _ 


acres 

acre-feet 

cubic  feet  per  second 
degrees  (angle) 
feet 

miles  (U.  S.  statute) 


4046*873 

1233.482 

0.02831685 

0.01745329 

0.3048 

1.609344 


METRIC  (SI) 

report  can  be  converted 

_ To  Obtain _ 

square  metres 

cubic  metres 

cubic  metres  per  second 

radians 

metres 

kilometres 
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JET  ENTRAINMENT  IN  PUMP ED-STO RAGE  RESERVOIRS 


1.  INTRODUCTION 

The  quality  of  water  extracted  from  a  reservoir  depends  strongly 
on  the  vertical  temperature  structure  existing  in  that  reservoir. 

Consider,  for  example,  the  interrelated  typical  water  quality  parameters, 
temperature  and  dissolved  oxygen.  As  a  result  of  surface  heating  by 
the  sun  and  vertical  mixing,  reservoirs  are  usually  thermally  stratified, 
with  the  temperature  decreasing  downwards  from  the  surface.  This  causes 
a  dynamically  stable  density  profile  which  reduces  vertical  mixing  by 
the  suppression  of  vertical  motions.  The  more  rapidly  that  ambient  density 
changes  with  depth,  the  more  that  these  motions  are  suppressed,  and  hence 
the  more  that  the  vertical  transfer  rates  are  suppressed.  This  can  cause 
the  rate  of  vertical  transfer  of  oxygen  to  be  so  low  that  oxygen  is 
depleted  faster  than  it  is  received  in  the  deeper  waters,  leading  to 
depressed  dissolved  oxygen  levels  there. 

This  pattern  can  be  significantly  different  in  pumped- storage 
reservoirs.  For  example,  if  the  inflow  of  pumpback  or  generation  is 
discharged  as  a  submerged  jet,  and  is  wanner  than  the  water  at  its  level 
of  injection,  buoyancy  forces  acting  on  the  jet  drive  it  towards  the 
surface.  The  combination  of  the  momentum  and  buoyancy  of  the  jet 
causes  it  to  entrain  ambient  reservoir  water.  Eventually,  this  rising  column 
will  reach  either  the  water  surface  or  a  level  of  equal  density,  at 
which  point  it  collapses  vertically  and  spreads  laterally.  Thus,  the 
process  has  redistributed  vertically  the  water  in  the  reservoir,  bringing 
denser,  deeper  water  towards  the  surface  and  altering,  both  directly 
and  indirectly,  the  dissolved  oxygen  content.  The  direct  change  is 
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the  redistribution  of  dissolved  oxygen  in  the  reservoir  by  this  vertical 
advection.  The  indirect  changes  are  due  to  the  altered  temperature 
structure  of  the  reservoir,  causing  changes  in  the  vertical  mixing 
coefficients  and  hence  changes  in  the  dissolved  oxygen  content.  Depending 
on  whether  water  is  withdrawn  from  near  the  surface  or  bottom  of  the 
reservoir,  these  changes  could  be  either  good  or  bad.  Thus,  the  tem¬ 
perature  and  dissolved  oxygen  as  well  as  other  water  quality  parameters 
of  water  withdrawn  from  a  pumped-storage  reservoir  will  generally  be 
different  from  those  withdrawn  when  the  reservoir  is  subject  only  to 
natural  temperature  variations. 

Prediction  of  the  change  in  vertical  temperature  distribution 
resulting  from  pumped-storage  operations  is  therefore  a  matter  of 
considerable  importance.  For  example,  it  may  be  necessary  to  predict  the 
performance  of  proposed  projects  to  assess  their  environmental  impact,  or 
to  aid  in  the  design  of  the  release  structures,  or  to  predict  the 
behavior  of  an  existing  project  under  different  operating  conditions. 

A  mathematical  model  of  general  applicability,  WESTEX,  has  been 
developed  by  the  U.  S.  Army  Corps  of  Engineers  Waterways  Experiment 
Station  (WES).  The  model  simulates  changes  in  temperature  stratification 
caused  by  pumped-storage  operations.  Vertical  mixing  is  simulated  by 
the  three  processes  of  entrainment,  mixing  of  inflow  and  entrained  flow, 
and  placement  of  the  mixture  at  the  appropriate  reservoir  level.  The 
model  requires  as  input  the  amount  of  entrainment  and  vertical  reservoir 
location  from  which  the  entrainment  occurred.  Entrainment  is  currently 
expressed  as  a  simple  multiple  of  the  incoming  flow,  where  the  multiple 
is  based  on  hydraulic  model  studies. 
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The  entrained  volume  depends  upon  the  momentum,  buoyancy,  and  volume 
fluxes  of  the  source,  the  discharge  structure  geometry,  and  the  ambient 
density  stratification.  Many  studies  have  been  made  of  entrainment  in 
different  fluid  flow  situations,  and  various  models  exist  for  the  pre¬ 
diction  of  entrainment.  The  present  study  was  commissioned  by  WES  to 
review  the  literature  pertaining  to  jet  entrainment  and  to  identify  math¬ 
ematical  descriptions  of  entrainment  which  offer  promise  in  improving 
mathematical  modeling  of  pumped- storage  reservoirs.  WES  will  investigate 
candidate  descriptions  in  a  subsequent  study.  The  results  of  the  review 
are  presented  in  this  report. 

The  report  is  divided  into  five  chapters.  In  the  following  chapter,  a 
discussion  of  typical  pumped-storage  projects  is  presented  to  illustrate 
the  types  of  flow  under  consideration.  In  Chapter  3,  the  important 
aspects  of  the  hydrodynamics  of  pumped-storage  reservoirs  are  discussed 
and  the  present  problem  defined.  Chapter  4  presents  the  methods 
recommended  for  entrainment  predictions  based  upon  the  literature  survey. 
The  flows  are  first  classified  by  a  primary  and  secondary  system  which 
is  explained  in  Section  4.1.  The  primary  classification  consists  of  four 
groups,  each  of  which  is  discussed  in  Sections  4.2  through  4.5.  The 
secondary  classification  also  has  four  parts,  each  of  which  is  discussed  in 
these  sections.  Each  of  these  16  types  of  flow  is  reviewed  separately 
and  recommendations  presented  for  the  prediction  of  entrainment  in  each 
case.  Finally,  in  Chapter  5  a  summary  of  the  study  and  conclusions  derived 
are  presented. 
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2.  TYPICAL  PUMPED- STORAGE  PROJECTS 


Pumped-storage  projects  exist  in  a  wide  variety  of  designs,  sizes, 
and  modes  of  operation,  causing  a  multitude  of  possible  inflow  character¬ 
istics.  The  type  of  system  considered  here  is  shown  schematically  in 
Figure  2.1;  it  consists  of  an  upper  and  a  lower  reservoir,  each  having  a 
free  surface.  The  reservoirs  are  connected  by  tunnels  or  pipes  within 
which  are  generating  and  pumping  units.  During  periods  of  peak  electric 
power  demand,  water  is  released  from  the  upper  reservoir  to  generate 
power  from  turbines.  When  power  demand  is  low,  for  example  at  night,  the 
flow  is  reversed  and  water  is  pumped  from  the  lower  to  upper  reservoir. 

Many  varieties  exist  within  this  basic  scheme.  Either  the  upper  or 
the  lower  reservoir  may  be  much  larger  than  the  other,  and  one  reservoir 
may  even  be  a  river.  The  intake/discharge  structures  can  have  various 
geometries,  and  may  have  separate  channels  for  inflow  and  outflow.  The 
discharge  may  be  near  either  the  bottom  or  the  surface  of  the  receiving 
water.  The  flow  into  the  reservoir,  which  may  be  either  to  the  lower 
reservoir  during  generation  or  to  the  upper  reservoir  during  purapback,  is 
of  interest  here.  Several  typical  pumped-storage  projects  are  outlined 
below  to  illustrate  the  diversity  of  flows  that  must  be  considered. 

2.1  Richard  B.  Russell  Lake 

This  project  is  described  in  Fontane  and  Bohan  (1974).  The  lake  im¬ 
pounds  the  Savannah  River  between  Clarks  Hill  Reservoir  and  the  Hartwell 
Dam  and  is  located  on  the  Georgia-South  Carolina  border.  The  Richard  B. 
Russell  Lake  forms  the  upper  reservoir  and  Clarks  Hill  Reservoir  the  lower 
when  the  system  is  used  in  the  pumped-storage  mode.  The  maximum  generation 
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Figure  2.1.  Schematic  diagram  of  pumped-storage  systems 


flow  from  Richard  B.  Russell  Lake  is  60,000  cfs*  and  the  maximum  pumpback 
flow  to  the  lake  is  20,000  cfs. 

2. 2  Marysville  Lake,  California 

This  proposed  project,  described  in  Fontane,  et  al.  (1977),  is  to  be 
situated  in  northern  California  on  the  Yuba  River.  The  upper  reservoir, 
Marysville  Lake,  would  be  a  multipurpose  lake.  It  would  be  much  larger 
than  the  lower  reservoir  with  a  storage  volume  varying  from  273,000  to 
916,000  acre-feet  and  a  surface  area  varying  from  3,180  to  6,640  acres. 


*  A  table  of  factors  for  converting  U.  S.  customary  units  of  measure¬ 
ment  to  metric  (SI)  units  is  presented  on  page  4. 


The  lower  reservoir  would  have  a  storage  volume  of  40,400  acre-feet  and 
an  area  of  947  acres.  Typically,  generation  and  pumpback  would  occur 
Monday  through  Saturday,  with  a  generation  flow  of  37,700  cfs  for  3.71 
hours  per  day,  and  a  pumpback  flow  of  8,200  cfs  for  9  hours  per  night. 

Initially  six  turbines  were  proposed  for  a  peak  generating  capacity  of 
900  MW. 

The  Waterways  Experiment  Station  physical  and  mathematical  model  studies 
were  concerned  with  definition  of  the  effects  of  pumpback  into  Marysville  Lake, 
whose  intake/discharge  structure  configuration  is  shown  in  Figure  2.2.  Out¬ 
flow  during  generation  could  be  through  any  of  the  five  selective  withdrawal 
intakes.  Each  intake  consisted  of  six  ports  50  feet  wide  by  10  feet  high  to 
accomodate  each  of  the  six  turbines.  Pumpback  flow  was  through  either  four  of 
the  selective  withdrawal  ports  or  a  28-foot  square  floodgate  near  the  bottom. 

Recently,  the  generating  capacity  of  this  system  was  increased  to 
2,250  MW  (Dortch  1978a),  although  most  of  the  project  features  remained 
unchanged.  The  previous  maximum  generation  flow  of  50,000  cfs  was  increased 
to  105,000  cfs  and  the  maximum  pumpback  flow  increased  from  9,000  to 
50,000  cfs.  The  storage  volumes  were  increased  somewhat,  and  five 
turbines  are  now  contemplated,  each  with  vertically  staggered  selective 
withdrawal  ports  25  feet  high  by  100  feet  wide. 

2.3  Dickey-Lincoln  School  Lakes,  Maine 

This  proposed  project  is  to  be  located  on  the  St.  John  River  in  north¬ 
ern  Maine  (Dortch,  et  al.  1976).  Again,  the  upper  reservoir,  Dickey  Lake, 
would  be  larger  than  the  lower  one,  extending  about  45  miles  upstream.  The 
maximum  storage  volume  of  this  lake  would  be  7,700,000  acre-feet,  with  a 
surface  area  of  86,000  acres.  The  lower  reservoir,  Lincoln  School  Lake, 
will  have  a  total  storage  of  about  85,000  acre-feet.  The  intakes  for 

10 


FIVE  SELECTIVE  WITHDRAWAL 
MTAKES.  EACH  INTAKE 


Figure  2.2.  Intake /discharge  and  reservoir  configurations  for 
Marysville  Lake  (900-MW  capacity) 
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generation  will  consist  of  a  multiple  penstock  intake  structure  with 
selector  gates  for  water  quality  control.  Flows  during  generation 
range  up  to  60,000  cfs.  Two  modes  of  operation  for  pumpback  are 
envisaged;  they  are  shown  in  Figure  2.3.  The  modes  are:  pumping 
over  the  top  of  the  selector  gates  (weirs)  and  pumping  under  the 
raised  gates. 

2.4  Allegheny  Reservoir,  Pennsylvania 

The  Allegheny  Reservoir-Kinzua  Dam  Project  is  located  on  the  Allegheny 
River,  Pennsylvania  (Fontane  and  Dortch  1978,  Dortch  1978b).  In  this 
case,  the  upper  reservoir  is  the  smaller  of  the  two;  it  has  a  capacity  of 
7,100  acre-feet  and  a  surface  area  of  110  acres.  The  lower  reservoir  is 
formed  by  the  impoundment  of  the  Allegheny  River  by  the  Kinzua  Dam.  Its 
capacity  is  about  572,000  acre-feet,  its  surface  area  12,080  acres, and 
its  depth  128  ft.  The  electrical  capacity  is  380  MW. 

2.5  Prattsville,  New  York 

This  proposed  project  would  be  about  40  miles  f  :  om  Albany,  New  York 
(Adams,  et  al.  1979).  The  upper  reservoir,  Schoharie  Reservoir,  would 
be  smaller  with  a  maximum  storage  of  30,000  acre-feet.  The  lower 
reservoir  would  have  a  maximum  storage  of  70,000  acre-feet.  The  generating 
capacity  would  be  1,000  MW,  the  maximum  generating  flow  rate  would  be  16,000 
cfs,  and  the  maximum  pumpback  flow  rate  would  be  10,668  cfs. 

2 . 6  Twin  Lakes,  Colorado 

These  lakes  are  located  at  an  elevation  of  about  9,000  feet  in 
central  Colorado;  together  they  constitute  the  lower  reservoir  in  this 
pumped-storage  project.  They  form  the  larger  of  the  reservoirs,  having 
a  combined  volume  of  170,000  acre-feet,  a  surface  area  of  about  3,000 
acres,  and  a  depth  of  nearly  100  feet.  The  upper  reservoir  has  a 
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volume  of  11,000  acre-feet*  Ultimately,  generation  and  pumpback  are 
contemplated  foi  Monday  through  Friday  with  a  generation  flow  of  7,200 
cfs  for  8  hours  per  day  and  a  pumpback  flow  of  6,180  cfs  for  9  hours 
per  night*  Two  units  are  proposed  for  a  capacity  of  200  MW.  For  a 
discussion  of  this  project,  see  King  and  Rhone  (1977). 
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3.  HYDRODYNAMICS  OF  PUMPED- STORAGE  RESERVOIRS 


3. 1  Introduction 

The  flows  and  mixing  induced  in  pumped-storage  reservoirs  involve 
unsteady,  three-dimensional,  density-stratified,  turbulent  processes 
for  which  no  general  predictive  method  exists*  These  processes  and 
studies  of  pumped-storage  reservoirs  are  briefly  reviewed  in  this  chap¬ 
ter  in  order  to  define  the  problem  under  consideration  and  present  the 
assumptions  made  in  solving  it. 

3.  ^  General  Features 

The  most  important  hydrodynamic  features  of  reservoirs  subject  to 
jetting  flows  are  illustrated  in  Figure  3.1.  The  jet  can  be  a  result 
of  pumpback  or  generation,  depending  on  whether  an  upper  or  a  lower 
reservoir  is  involved. 

Reservoirs  subjected  to  natural  temperature  variations  usually  become 
density  stratified  during  spring  and  summer,  as  shown  in  Figure  3.1a, 

Maximum  stratifications  occur  during  summer.  Minimum  stratifications  occur 
during  winter  when  the  reservoir  may  be  well-mixed.  Usually,  there  are  only 
small  lateral  currents,  and  the  shape  of  the  density  profile  will  be  similar 
at  different  locations  within  the  reservoir. 

The  flow  field  shortly  after  the  start  of  generation  or  pumpback 
for  the  case  of  a  positively  buoyant  jet  is  shown  in  Figure  3.1b. 

Gravity  forces  resulting  from  the  buoyancy  cause  the  trajectory  of  the 
jet  to  curve  upwards.  As  the  jet  ascends,  it  entrains  fluid  from  the 
surrounding  reservoir,  and  the  volume  flow  rate  within  the  jet  increases 
with  distance  from  its  origin.  If  the  reservoir  is  density  stratified. 
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Initial  conditions, 

typical  density  stratification 


b.  Buoyant  Jotting  billow 


c.  Shortly  altar  billow  caaaaa 


cL  Long  tbna  altar  billow  caaaaa 


Figure  3.1  Important  hydrodynamic  features  of  pumped- storage 
reservoirs  subject  to  jetting  inflows 


1 

1 

] 
\ 
i 

the  jet  can  eventually  reach  a  level  at  which  its  density  is  equal  to 
that  of  the  background  water.  At  this  point,  the  jet  collapses  vertically 
and  begins  to  spread  horizontally.  If  the  ambient  stratification  is  not 
strong  enough  to  keep  the  jet  submerged,  it  will  reach  the  surface.  In 
this  situation  the  deeper  waters  are  "pumped"  upwards,  causing  vertical 
mixing  within  the  reservoir. 

Shortly  after  flow  ceases  (Figure  3.1c)  the  mixture  of  inflow  and 
ambient  water  continues  to  spread  horizontally  as  a  density  current, 
during  which  mixing  may  continue  within  the  layer.  Eventually  the  layer 
will  reach  the  reservoir  boundaries,  and  reflections  may  occur.  Finally, 
however,  horizontal  motions  should  essentially  cease,  and  the  density 

! 

profile  again  become  similar  at  different  locations  in  the  reservoir.  The 

t 

density  profile  at  this  time  will  be  different  from  the  original  one 

i 

(Figure  3. Id),  the  effect  of  the  process  being  to  make  the  density  profile 
more  uniform  with  depth. 

The  above  discussion  illustrates  only  the  essential  features  of 
importance  here.  Other  features  of  secondary  importance,  for  example, 
the  generation  of  internal  waves  in  the  reservoir,  will  not  be  discussed. 

i 

Also  the  horizontal  density  current  may  not  reach  the  reservoir  boundaries  j 

before  withdrawal  from  the  reservoir  begins,  nor  even  before  the  next 
cycle  begins.  Therefore,  steady  state  conditions  may  never  be  reached 
in  the  reservoir  during  pumped-s torage  operations. 
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The  present  study  is  concerned  only  with  the  entrainment  during 
inflow,  the  situation  shown  in  Figure  3.1b.  A  detailed  definition  sketch 
of  the  problem  is  shown  in  Figure  3.2  for  a  positively  buoyant  discharge. 

Conservation  of  volume  for  an  incompressible  fluid  yields: 

Qc  "  Qo  +  % 

where  Qc*  is  the  total  volume  flow  rate  in  the  density  current,  Qq  is  the 
volume  flow  rate  of  the  inflow,  and  Qg  is  the  total  volume  flow  rate  of 
entrainment  of  ambient  water  by  the  jet.  An  overall  entrainment  coeffi¬ 
cient,  E,  is  now  defined  by: 

Qe  -  EQC  (3.2) 

The  purpose  of  this  study  is  to  provide  a  method  of  predicting  E  and  the 
variation  of  entrainment  with  depth. 

3.3  Hydraulic  Model  Studies 

Only  a  few  hydraulic  model  studies  of  pumped-storage  reservoirs  of 
particular  geometries  or  of  general  applicability  have  been  performed  to 
investigate  the  effects  of  pumpback  and  generation.  Some  of  these  studies 
are  summarized  below. 

Fontane  and  Bohan  (1974)  conducted  studies  on  the  effects  of  pumpback 
flow  during  pumped-storage  operation  on  the  Richard  B.  Russell  Lake.  They 
found  an  entrainment  coefficient,  E,  of  1.5,  where  67  percent  of  the  en¬ 
trained  flow  was  from  the  hypolimnion  and  33  percent  from  the  epilimnion. 
The  complexity  of  the  entrainment  process  was  illustrated  by  the  fact  that 
following  the  changeover  from  generation  to  pumpback,  the  level  of 


*  For  convenience,  symbols  are  listed  and  defined  in  the  Notation 
(Appendix  A). 
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Figure  3.2.  Definition  sketch  of  inflow 
jet  entrainment 

entrainment  was  limited  to  a  layer  near  the  intake  level  and  not  to  the 
bottom  reservoir  levels. 

Dortch,  et  al.  (1976)  describe  the  results  of  studies  on  selective 
withdrawal  and  inflow  in  Dickey-Lincoln  School  Lakes.  The  inflow  tests 
were  conducted  for  flows  over  and  under  the  gates  (see  Figure  2.3). 
Entrainment  coefficients,  E,  were  found  to  be  0.3  for  pumping  over  the 
gates  for  both  negatively  and  positively  buoyant  jets  and  0.6  for  pumping 
under  the  gates  with  a  positively  buoyant  jet.  Stratification  in  Lincoln 
School  Lake  was  completely  destroyed  within  two  simulated  weeks  of  opera¬ 
tion.  After  this  period  of  time  the  larger  reservoir  (Dickey  Lake)  was 
still  strongly  stratified,  although  some  vertical  mixing  had  occurred. 
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The  depth  from  which  the  entrained  flow  originated  depended  on  the 
flow  situation.  For  a  negatively  buoyant  inflow  over  the  weir,  the 
entrained  flow  came  primarily  from  above  the  weir  elevation.  For  positive 
buoyancy,  the  entrainment  was  from  all  depths  below  the  weir  elevation, 
with  the  maximum  contribution  at  20  feet  below  the  weir  crest.  For 
positively  buoyant  flow  from  below  the  gates,  entrainment  was  primarily 
from  the  deeper  waters.  In  all  of  these  cases,  the  layers  contributing 
to  the  entrained  flow  depend  upon  the  particular  density  stratification 
profile  present. 

Fontane,  et  al.  (1977)  report  similar  tests  performed  for  Marysville 
Lake,  California.  The  discharge  ports  in  this  situation  are  shown  in 
Figure  2.2.  For  pumpback  from  the  upper  ports,  E  was  found  to  be  about 
0.7.  The  entrained  flow  came  primarily  from  the  levels  at  and  immediately 
above  the  port  level.  For  pumpback  from  the  floodgate,*  E  increased  to 
between  2.5  and  3.0,  with  the  entrained  flow  mainly  from  the  deeper  waters. 
This  increased  entrainment  is  due  to  the  positive  buoyancy  of  the  inflow. 

Tests  were  also  conducted  for  this  lake  when  the  electrical  capacity 
of  the  project  was  increased  from  900  to  2250  MW,  with  a  corresponding 
increase  in  flow  rates  and  port  dimensions.  The  results,  presented  in 
Dortch  (1978a),  suggested  an  increase  in  E  to  1.0  for  discharge  through  an 
upper  port.  The  entrained  flow  was  again  observed  to  come  mainly  from 
the  levels  of  the  port. 

Tests  on  the  proposed  pumped  storage  plant  at  Prattsville  were 
reported  by  Brocard  and  Nyquist  (1978).  Values  of  entrainment  coeffi¬ 
cients,  E,  were  evaluated  indirectly  from  these  results  by  Adams,  et  al. 
(1979),  and  a  value  of  E  =  1  was  chosen. 
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Few  studies  have  been  performed  from  which  results  of  general 
applicability  can  be  derived.  Zuluaga- Angel ,  et  al.  (1972)  and 
Darden,  et  al.  (1975)  report  on  the  circulation  created  by  a  turbulent 
slot  jet  discharging  into  a  linearly  stratified  rectangular  reservoir. 

They  find  complex  circulation  patterns  to  be  induced,  with  no  true 
steady  state  flow  existing.  Entrainment  increased  considerably  with  the 
buoyancy  of  the  inflow. 

3. 4  Mathematical  Models 

Several  mathematical  models  have  been  developed  to  predict  temperature 
stratifications  due  to  pumped-storage  reservoir  operations.  All  of  these 
models  are  one-dimensional  in  their  treatment  of  the  reservoir,  that  is, 
the  reservoir  consists  of  a  series  of  horizontal,  well  mixed  layers;  the 
treatment  of  the  inflow  jet  may  be  three-dimensional.  These  models  are 
summarized  below. 

The  M. I.T.  model  (Adams,  et  al.,  1979)  predicts  the  vertical 
distribution  of  reservoir  temperature.  The  model  accounts  for  heat 
fluxes  to  each  horizontal  element  due  to  solar  radiation,  surface  heat 
transfer,  inflows  and  outflows,  turbulent  diffusion,  surface  wind  mixing, 
and  convective  overturning.  The  entrainment  due  to  pumped-storage 
operations  is  assumed  in  this  case  to  come  from  the  layers  above  the 
intake  and  to  be  uniformly  distributed  over  these  layers.  By  comparing 
the  results  of  mathematical  simulations  with  the  results  of  physical 
model  tests,  a  value  of  E  =  1  was  assumed.  The  principal  finding  of  the 
results  of  simulations  of  the  proposed  reservoir  was  a  mixing  of  the 
upper  layers  with  a  deepening  of  the  epilimnion.  The  simulation  results 
were  found  to  be  insensitive  to  the  value  of  E  in  the  range  0.5  <  E  <2. 

The  model  was  also  used  to  simulate  temperatures  in  the  operating  Blenheim- 


Gilboa  Station;  reasonable  agreements  between  predicted  and  measured 
results  were  obtained,  with  most  discrepancies  occurring  in  spring, 
when  surface  temperatures  were  overestimated. 

A  similar  model  has  been  developed  by  the  U.S.  Army  Corps  of  Engineers 
(Dortch,  et  al .  1976).  The  entrainment  coefficient,  E,  is  based  on 
the  results  of  a  physical  model  study  for  each  project,  as  discussed  in 
Section  3.3.  The  entrained  flow  is  assumed  to  come  entirely  from  one 
layer  of  the  reservoir,  the  layer  being  chosen  from  the  results  of  the 
physical  model  study.  The  average  density  of  the  mixture  is 
then  computed  and  the  flow  placed  at  its  level  of  neutral  buoyancy. 
Subsequent  modifications  to  the  model  are  described  in  Dortch  (1978a), 
where  computation  of  entrained  flow  from  several  layers  is  allowed, 
improving  the  degree  of  realism  of  the  model.  Applications  of  the 
model  to  the  Allegheny  Reservoir  (Dortch  1978b)  produced  good  simula¬ 
tions  of  the  observed  temperature  profiles.  Further  refinement  of  the 
entrainment  description  (Fontane  and  Dortch  1978)  produced  even  closer 
simulations  to  the  observed  temperatures. 

A  model  was  developed  by  Ditmars  (1970)  to  study  mixing  in  stratified 
reservoirs.  In  the  situation  considered,  mixing  was  intentionally  induced 
by  extracting  warm  water  from  near  the  surface  and  reinjecting  it  as  a 
positively  buoyant  jet  near  the  bottom.  Thus,  the  model  was  not  intended 
specifically  for  pumped-storage  applications,  although  it  can  be  easily 
adapted  to  this  case.  A  three-dimensional  integral  analysis  of  the  jet 
is  used,  and  entrainment  into  the  rising  jet  is  allowed  continuously 
from  each  layer.  Ditmars  presents  general  results  from  which  estimates 
can  be  made  of  the  time  required  for  complete  mixing  of  a  linearly 
stratified  reservoir.  The  discharge  of  the  warm  upper  layers  of  water 


as  a  buoyant  jet  into  the  cold  lower  layers  maximizes  entrainment  and 
hence  maximizes  the  vertical  mixing  rate.  Therefore,  the  results  can 
be  used  to  provide  a  lower  bound  for  the  time  required  for  complete 
mixing  of  a  pumped-storage  reservoir.  This  time  is  proportional  to  the 
ratio  of  the  reservoir  volume  between  the  intake  and  discharge  to  the 
pumping  rate. 

A  similar  model  is  presented  by  Henderson-Sellers  (1978). 

3.5  Conclusions 

The  amount  of  mixing  induced  in  pumped-storage  reservoirs  can  vary 
widely,  ranging  from  complete  mixing  induced  in  time  scales  on  the  order 
of  a  week  to  small  changes  in  stratification  over  months  of  operation. 

In  general,  the  smaller  the  reservoir  relative  to  the  daily  amount  of 
water  pumped,  the  shorter  will  be  the  mixing  time.  The  mixing  time 
also  depends  on  the  amount  of  entrainment;  discharge  structure  designs 
which  increase  the  amount  of  entrainment  will  decrease  the  mixing  time. 
Thus,  the  probable  reason  for  the  apparent  insensitivity  to  E  shown  by 
some  of  the  mathematical  simulations  discussed  above  is  the  large  size 
of  the  reservoirs  considered.  For  smaller  reservoirs  the  amount  of 
entrainment  can  have  a  critical  effect  on  the  density  stratification. 

The  process  of  entrainment  is  a  complex  one;  the  location  of  the 
layers  from  which  the  entrainment  comes  is  a  function  of  the  shape  and 
strength  of  the  ambient  density  stratification,  as  well  as  the  buoyancy 
and  momentum  fluxes  of  the  discharge  and  its  location.  As  pumpback  or 
generation  occurs  for  only  a  few  hours,  entrainment  is  assumed  constant 
during  this  time,  and  the  effects  of  reservoir  geometry  are  assumed  negli 
gible.  For  continuous  discharge  these  assumptions  would  not  be  valid. 
Methods  for  predicting  entrainment  amounts  and  sources  using  these  assump 
tions  are  discussed  in  the  following  chapter. 
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4.  PREDICTION  OF  JET  ENTRAINMENT 


4 . 1  Introduction 

The  many  parameters  affecting  entrainment  in  pumped-storage  reservoirs 
result  in  a  complex  problem  of  which  many  studies  have  been  made.  The 
relevant  studies  are  reviewed  in  this  chapter  and  potential  techniques 
for  entrainment  prediction  are  presented. 

Entrainment  occurs  as  a  result  of  turbulent  shear  flows.  As  discussed 
by  Roshko  (1976),  entrainment  is  the  incorporation  of  nonturbulent ,  usually 
irrotational  fluid  into  the  turbulent  region  or,  conversely,  the  diffusion 
of  the  turbulent  region  into  the  ambient  flow.  Entrainment  appears  to  be 
due  to  the  engulfing  action  of  large  turbulent  eddies  near  the  jet 
boundaries.  Free-stream  fluid  is  drawn  in  between  vortices  and  ingested 
into  the  jet  where  it  is  made  turbulent  and  digested  by  the  action  of 
smaller  eddies.  The  rate  of  the  entrainment,  therefore,  depends  on  the 
nature  of  these  large  eddies.  These  depend  on  the  source  geometry,  the 
momentum  and  buoyancy  of  the  discharge,  ambient  stratification  and 
proximity  to  local  boundaries. 

As  entrainment  is  local,  the  overall  reservoir  size  and  geometry 
should  not  affect  the  entrainment  rate  and  their  effects  will  be  assumed 
negligible.  The  size  and  geometry  will,  however,  affect  overall  circula¬ 
tion  patterns  in  the  reservoir,  and  possibly  the  source  of  the  entrained 
flow,  For  the  short  times  over  which  flows  occur,  entrainment  into  the  jet 
will  be  assumed  to  be  a  steady  process;  this  assumption  would  be  invalid  for 
longer  times.  As  the  entrainment  depends  more  on  the  source  buoyancy  and 
momentum  fluxes  than  on  the  volume  flux,  the  effect  of  source  volume  flux 


can  usually  be  neglected  within  the  accuracy  of  engineering  estimations. 

Some  of  the  predictive  techniques  presented  in  this  chapter  neglect  the 
effect  of  source  volume  flux  on  entrainment,  but  some  do  not. 

Even  with  these  simplifying  assumptions,  no  general  technique  exists  for 
prediction  in  all  possible  solutions.  Consider  the  types  of  flows 
considered  in  Chapter  2,  for  example.  The  jet  can  be  near  the  surface  or 
the  bottom  (Figure  2.3).  It  can  be  from  a  structure  approximating  a  slot 
or  square  orifice  (Figure  2.2).  Depending  on  where  in  the  other  reservoir 
the  jet  came, the  jet  could  be  heavier  or  lighter  than  or  of  equal  density 
to  the  receiving  water.  And  finally,  depending  on  the  reservoir  size  and 
duration  of  pumpback  or  generation,  the  receiving  waters  can  be  density 
stratified  or  well  mixed.  Therefore,  to  account  for  all  the  types  of  flow 
which  could  be  conceived,  it  is  necessary  to  devise  a  scheme  for  their 
class if icat ion . 

The  flows  are  classified  according  to  a  primary  and  secondary  system.  The 
primary  classification  scheme  depends  on  the  source  buoyancy  and  boundary  prox¬ 
imity;  it  is  shown  in  Figure  4,1,  in  which  four  classes  exist.  The  secondary 
classification  scheme  depends  on  the  ambient  stratification  and  source  geometry. 
The  ambient  water  is  either  stratified  or  uniform,  and  the  discharge  approximates 
a  slot  or  round  orifice.  Thus,  each  of  the  four  flow  schemes  shown  in  Figure 
4.1  has  four  subgroups  for  a  total  of  16  flow  types.  Each  oi  these  16  flow- 
types  is  discussed  in  this  chapter. 

The  chapter  is  divided  according  to  the  flow  classification  scheme.  Each 
primary  group  is  discussed  in  Sections  4.2  through  4.5.  Each  of  these  sections 
is  further  divided  into  four  parts  where  each  flow  type  is  discussed.  The 
parts  consist  of  a  brief  review  of  pertinent  work,  followed  by  a  recommended 
predictive  technique.  The  techniques  exist  in  either  analytical  form  as 


25 


b.  Nonbuoyant  bniMTMd  Jot 


c.  Buoyant  boundary  Jot 


s- 


d.  Nonbuoyant  boundary  Jot 

Figure  4.1.  Primary  flow  classification  scheme 
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equations,  as  graphs,  or  as  mathematical  models. 

In  many  of  the  studies  discussed  below,  entrainment  is  described  indi¬ 
rectly  by  dilution.  Consider  a  mixture  consisting  of  a  volume  of  source 
flow  mixed  with  a  volume  of  receiving  water.  The  dilution,  S,  is  defined 


V  +  V 
o  r 


(4.1) 


Similarly,  dilutions  can  be  related  to  volume  fluxes  Qq  and  ;  this 
will  be  done  in  the  following  sections  if  solutions  are  given  in  terms 
of  dilution. 

4.2  Buoyant  Immersed  Jet 

This  situation  is  case  a.  Figure  4.1.  It  occurs  frequently  in 
practice,  as  heated  discharges  from  power  plants  into  lakes  and  coastal 
waters,  and  domestic  sewage  discharges  into  coastal  waters.  Consequently, 
many  studies  of  the  phenomena  have  been  made  and  predictive  techniques 
for  the  various  flows  exist.  These  flows  consist  of  horizontal  discharges 
from  slot  or  round  orifices  into  a  stratified  or  unstratified  fluid. 

They  are  reviewed  in  this  section. 

A  definition  sketch  showing  the  details  of  the  problem  considered  is 
shown  in  Figure  4.2,  The  discharge  is  horizontal  and  the  ambient  fluid 
may  or  may  not  be  density  stratified.  The  density  of  the  discharged 
water,  p^,  is  greater  than  or  less  than  the  density  of  the  ambient  fluid, 
at  the  level  of  the  discharge,  causing  buoyancy  forces  to  act  on  the 
jet.  The  particular  situation  shown  in  Figure  4.2  is  for  the  case  of  a 
positively  buoyant  (p^  <  Pq)  discharge  into  a  stratified  fluid.  The 
ambient  fluid  density,  pfl,  is  a  function  of  depth,  y.  For  a  round  orifice 
the  diameter  is  D,  and  for  a  slot  the  width  is  B.  The  governing  source 
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parameters  can  be  defined  as  the  source  volume,  momentum  and  buoyancy 

fluxes:  Q  ,  M  »  and  B  for  a  round  discharge;  and  q  ,  m  ,  and  b  per  unit 
o  o  o  ^o  o  o 

width  for  a  slot  discharge,  where: 


IT  _2 

TlID 

4 

(4.2) 

UQ0 

(4.3) 

(4.4) 

uB 

(4.5) 

(4.6) 

(4.7) 

VP1 

p0  8 

0 

(4.8) 

u  is  the  mean  initial  jet  velocity,  and  g  the  acceleration  due  to  gravity. 

The  general  features  of  the  flow  are  also  shown  in  Figure  4.2. 

These  features  are  discussed  in  many  other  papers,  for  example  Koh  and 
Brooks  (1975).  Near  the  orifice  is  a  zone  of  flow  establishment  in 
which  the  velocity  profile  across  the  jet  width  becomes  established. 

This  zone  extends  for  about  6. 2D  for  round  jets  and  5.2B  for  slot  jets 
(Albertson,  et  al.  1950).  During  the  zone  of  flow  establishment  and 
beyond,  the  jet  entrains  ambient  fluid.  Near  the  source  the  entrain¬ 
ment  is  primarily  due  to  the  source  momentum  flux;  far  from  the  source 
the  buoyancy  forces  and  hence  the  source  buoyancy  flux  dominate  the 
entrainment.  Eventually  the  rising  fluid  column  either  reaches  the 
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water  surface  or  a  level  at  which  its  density  is  equal  to  the  ambient 
density.  (If  the  ambient  fluid  is  unstratified,  the  jet  will  always 
reach  the  surface.)  Having  reached  its  level  of  neutral  buoyancy,  the  flow 
field  collapses  vertically  and  begins  to  spread  horizontally.  The 
thickness  of  the  spreading  layer,  h,  can  be  estimated  from  the  results 
for  a  two-dimensional  line  source  of  buoyancy  flux  in  an  unstratified 
ambient  fluid  as  (Roberts  1979) : 
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The  functional  form  of  Eq.  4.11  is  given  in  Figure  4.3.  For  values  of 
discharge  Froude  number  less  than  8,  a  mathematical  model  discussed  in 
Section  4,2.2  should  be  used. 

In  order  to  use  Figure  4.3  to  estimate  E,  E  and  must  be  related. 
It  can  be  shown  (Brooks  1972,  Fan  and  Brooks  1969)  that  the  centerline 
dilution  is  given  by: 


m 


1.15 


Eqs.  3.1,  3.2,  and  4.12  are  combined  to  give 
S 

-  1. 


1.15 


(4.12) 


(4.13) 


Hence,  Figure  4.3  in  conjunction  with  Eqs.  4.12  and  4.9  can  be  used  to 
estimate  E.  The  total  height  of  rise  is  H  (see  Figure  4.2). 

This  technique  can  also  be  used  where  the  receiving  water  is  of 
uniform  density  up  to  the  epilimnion,  where  the  density  changes  suddenly 
with  a  magnitude  large  enough  to  trap  the  rising  plume  beneath  it.  In 
this  case  the  rise  height,  yfc ,  is  the  water  depth  below  the  epilimnion. 

The  layers  from  which  entrainment  occurs  can  be  estimated  by  computing 
E  for  various  heights  y  and  applying  Eq.  3.2. 

4.2.2  Round  discharge,  stratified  ambient 

This  situation  could  occur  in  the  Marysville  project  (Chapter  2) 
if  a  buoyant  inflow  was  discharged  through  the  28-foot  square  floodgate 
near  the  reservoir  bottom  (see  Figure  2.2).  The  square  floodgate  can  be 
approximated  by  a  round  port  of  equal  area,  and  the  techniques  of  this 
section  applied.  There  are  no  general  models  for  discharges  other  than 


those  from  round  or  slot  orifices. 


Dilution  and  entrainment  can  be  estimated  from  the  integral  solutions 
of  Fan  and  Brooks  (1969).  Their  report  contains  the  details  of  the 
solution  method  which  will  not  be  repeated  here.  Fan  and  Brooks  show 
that  for  linear  stratifications,  the  dilution  at  the  terminal  rise  height, 
St ,  can  be  expressed  as: 


S 


t 


(4.14) 


where  u  =  f (m  ,  u  ) 
t  o  o 


(4.15) 


m  and  jj  are  dimensionless  momentum  and  volume  flux  parameters  given 
o  o 

by  Brooks  (1972)  as: 


and 


m  =  0.324  F2T  1 
o 


„  -2.38  F1/4I'5/8 

o 


T  is  a  stratification  parameter  defined  as 


0.87  (p  -pj 

O  L 


D  (-dp  /dy) 


(4.16) 

(4.17) 

(4.18) 


The  functional  form  of  Eq.  4.15  is  shown  in  Figure  4.4.  This  figure  can 
be  used  in  conjunction  with  Eqs.  4.13  and  4.14  (with  =  S^_)  to  estimate 
E.  Again,  corrections  should  be  made  for  the  blocking  effect  (Eq.  4.9). 

The  terminal  rise  height,  y  ,  is  also  given  by  Fan  and  Brooks.  They 
show  that  y  is  given  by: 


1/43/8 


5*-  =  i-37  V  T 


(4.19) 


where 


r 

't 


is  a  dimensional  height  of  rise  and 

=  f(m  ,u  ) 
o  o 


(4.20) 
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Figure  4,4. 


Terminal  volume  flux  parameter  pt  for  horizontal  round 
buoyant  jets  (from  Fan  and  Brooks  1969) 
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The  functional  form  of  Eq.  4.20  is  shown  in  Figure  4.5. 

To  estimate  the  layers  from  which  entrainment  comes,  and  to  predict 
overall  entrainment  when  the  stratification  is  nonlinear,  it  is  best  to 
use  the  original  computer  programs  from  which  these  results  were  generated. 
The  model  is  described  in  Ditmars  (1969).  The  levels  from  which  entrainment 
comes  can  be  obtained  by  inspection  of  the  output,  which  gives  dilution 
as  a  function  of  depth,  or  by  slight  modification  of  the  code  to  give  the 
result  directly. 

4.2.3  Slot  discharge,  unstratified  ambient 

This  problem  has  been  analyzed  by  similar  integral  techniques  to 
those  discussed  in  Section  4.2.2.  Brooks  (1972)  shows  that  the  centerline 
dilution,  Sq,  can  be  expressed  as 

So  =  f  (b  ’  F)’  (A-21) 

where  F  is  a  Froude  number  defined  by 


F  =  —= —  (4.22) 

Vs'B 

and  u  is  the  velocity  of  efflux  at  the  source.  The  functional  form  of 
Eq.  4.21  is  shown  in  Figure  4.6.  For  a  slot  jet,  the  centerline  dilution 
is  equal  to  the  volume  flux  ratio,  or 


S 

o 


(4.23) 


Combining  Eqs.  3.1,  3.2,  and  4,23  gives 


E  =  Sq  -  1  (4.24) 

Hence,  Figure  4.6  can  be  used  in  conjunction  with  Eq.  4.24  to  estimate  E 
after  correction  is  made  for  blocking  (Eq.  4.9).  In  this  unstratified  case 
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Figure  4.5.  Terminal  rise  height  f,  t  for  horizontal  round  buoyant  jets 

(from  Fan  and  Brooks  1969) 


the  terminal  rise  height,  y^,  is  equal  to  the  water  depth,  d. 

This  technique  can  also  be  used  where  the  receiving  water  is  of 
uniform  density  up  to  the  epilimnion,  where  the  density  changes  suddenly 
and  with  a  magnitude  large  enough  to  trap  the  rising  plume  beneath  it. 

In  this  case  the  rise  height,  y  ,  is  equal  to  the  water  depth  beneath 
the  epilimnion. 

The  layers  from  which  entrainment  comes  can  be  estimated  by 
computing  E  for  various  heights  y  and  applying  Eq .  3.2. 

4.2.4  Slot  discharge,  stratified  ambient 

A  possible  occurrence  of  this  situation  would  be  the  Marysville 
project  if  a  buoyant  inflow  were  discharged  via  one  of  the  lower 
selective  withdrawal  intakes  (Figure  2.2).  Each  of  these  intakes  consists 
of  six  ports,  each  50  feet  long  by  10  feet  high  (for  the  900-MW  capacity). 
Hence,  if  the  inflow  was  warmer  than  the  local  ambient  water,  the  dis¬ 
charge  would  approximate  the  two-dimensional  buoyant  slot  jet  into  a 
stratified  fluid,  the  problem  discussed  in  this  section. 

The  situation  has  been  analyzed  by  the  integral  techniques  previously 
discussed.  Brooks  (1972)  shows  that  for  a  linear  stratification,  the 
centerline  dilution  at  the  terminal  rise  height,  Sq,  can  be  expressed  as: 

s  25) 

O  1  t  0 

u  =  f (m  , u  )  (4.26) 

t  o  o 

where  mQ  and  are  dimensionless  momentum  and  volume  flux  parameters 
given  by 

m  =  0.5  F2  T_1  (4.27) 

o 
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and 


(4.28) 


2/3  - 

P  =  1.85  F  '  T 
o 


F  is  defined  by  Eq.  4.22,  and 


T  = 


Po'Pl 


B(~dpa/dy) 


(4.29) 


The  functional  form  of  Eq.  4.26  is  shown  in  Figure  4.7,  which  can  be 
used  with  Eq.  4.24  to  estimate  E.  Again,  correction  for  blocking  should 
be  made  (Eq.  4.9) . 

The  terminal  rise  height,  yt,  can  also  be  estimated.  Brooks  (1972) 
shows  that 


~  =  0.96C  F1/3  T1/2 
j5  t 


(4.30) 


where  £  =  f (m  ,y  ) 

t  o  o 


(4.31) 


The  functional  form  of  Eq.  4.31  is  shown  in  Figure  4.8. 

To  estimate  the  layers  from  which  the  entrained  water  comes  and  to 
predict  gross  entrainment  when  the  stratif ication  is  nonlinear,  it  is 
best  to  use  the  original  computer  programs  from  which  these  results  were 
generated.  The  programs  are  given  in  Sotil  (1971).  The  levels  from 
which  entrainment  comes  can  be  obtained  by  inspection  of  the  output, 
which  gives  dilution  as  a  function  of  rise  height,  or  by  slight  modifica¬ 
tion  of  the  code  to  give  the  result  directly. 

Other  relevant  experiments  were  reported  in  Zuluaga-Angel,  et  al. 
(1972)  and  Darden,  et  al.  (1975).  They  measured  dilution  and  observed 
the  flow  field  caused  by  discharge  of  buoyant  and  nonbuoyant  slot  jets 
into  linearly  stratified  model  prismatic  reservoirs.  The  results  for 
nonbuoyant  discharges  are  discussed  in  Section  4.3.4.  For  buoyant 
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Figure  4.8.  Terminal  height  of  rise  for  horizontal  slot  buoyant  jets 
in  linearly  density  stratified  environments  (from  Fan  and  Rrooks 


discharges,  they  found  the  entrainment  to  be  much  greater  than  for 
nonbuoyant  discharges.  Although  no  method  is  given  for  predicting 
entrainment  of  a  buoyant  jet,  the  results  provide  useful  informa¬ 
tion  on  the  resulting  complex  reservoir  flow  patterns,  and  hence 
entrainment  layers. 

4 . 3  Nonbuoyant  Immersed  Jet 

This  situation  is  case  b.  Figure  4.1.  The  density  of  the  jet  is 
equal  to  the  ambient  density  at  its  level  of  injection,  so  no  buoyancy 
forces  exist.  For  the  purpose  of  engineering  calculations,  the  entrain¬ 
ment  can  be  assumed  to  be  caused  entirely  by  the  source  momentum  flux. 
This  assumption  leads  to  simple  analytical  solutions  in  the  unstratified 
case,  although  the  stratified  case  is  more  complex. 

4.3.1  Round  discharge,  unstratified  ambient 

The  volume  flux  in  the  jet,  Q^,  is  assumed  to  depend  only  on  the 
initial  momentum  flux,  M^,  and  the  distance,  x,  from  the  source,  thus 

Qc  =  f<M0»x>  (4.32) 

A  dimensional  analysis  of  Eq.  4.32  yields 


where 

round 


Cj  is  a  constant  which  must  be  measured  experimentally, 
port : 


(4.33) 


For  a 


u  » 


(4.34) 


and  a  combination  of  Eqs.  4.3,  4.33,  and  4.34  yields 
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(4.35) 


Qo 


X 

D 


Thus,  the  volume  flux  in  the  jet  increases  linearly  with  distance  from 
the  source.  The  experimental  work  of  Albertson,  et  al.  (1950)  yields 
-  0.284,  so  Eq.  4.35  becomes 


^  =  °'32  D 


(4.36) 


Substitution  of  Eqs.  3.1  and  3.2  into  Eq.  4.36  yields 


E  =  0.32  ^  -  1 


(4.37) 


It  is  important  to  note  that  in  this  case,  E  is  not  equal  to  a  finite 
value  but  increases  with  x,  and  therefore  time.  To  estimate  the  total 
volume  of  entrainment  for  the  duration  of  the  pumped-storage  inflow 
cycle,  Eq.  4.37  must  be  integrated  with  respect  to  time.  To  do  this, 
the  length  of  the  jet  as  a  function  of  time  is  required;  this 
can  be  obtained  from  the  results  of  Albertson,  et  al.  (1950). 

As  the  density  of  the  ambient  water  is  uniform,  the  density  of 
the  mixture  is  indifferent  to  the  location  of  the  entrained  flow. 

4.3.2  Round  discharge,  stratified  ambient 

The  effect  of  stratification  is  to  suppress  vertical  turbulent 
fluctuations,  leading  to  eventual  collapse  of  the  jet  and  suppression 
of  entrainment.  The  problem  is  analogous  to  that  of  wake  collapse 
behind  a  submarine  in  a  stratified  fluid,  reviewed  by  Lin  and  Pao  (1979). 

E  is  finite  and  should  depend  on  the  strength  of  the  stratification. 
No  experimental  data  were  found  from  which  E  would  be  predicted.  From 
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the  experimental  data  on  pumped-storage  reservoirs  discussed  in  Section 
3.3,  E  would  not  be  expected  to  differ  significantly  from  unity  for  con- 
ditions  typical  of  pumped-storage  projects. 


Because  of  the  suppression  of  vertical  motions  by  the  stable  strati¬ 
fication,  most  of  the  entrained  flow  should  come  from  the  same  level  as 
the  discharge  port. 

4.3.3  Slot  discharge,  unstratified  ambient 

The  volume  flux  is  again  assumed  to  depend  only  on  the  source 

momentum  flux,  m  : 

o 

q  =  f (m  ,x)  (4.38) 

c  o 

where  is  the  volume  flux  per  unit  width.  A  dimensional  analysis 
of  Eq.  4.38  yields: 


Again,  E  does  not  tend  to  a  finite  value  but  increases  continuously 
with  x,  and  therefore  time.  To  estimate  the  total  volume  entrained 
during  the  pumped-storage  inflow  cycle,  Eq.  4.41  must  be  integrated  with 
respect  to  time.  To  do  this,  the  length  of  the  inflow  jet  as  a  function 


of  time  is  needed;  this  can  be  obtained  from  the  experiments  of  Albertson, 
et  al.  (1950). 

As  the  density  of  the  inflow  and  ambient  water  is  uniform,  the 
density  of  the  jet  is  indifferent  to  the  location  of  the  entrained  flow. 
4.3.4  Slot  discharge,  stratified  ambient 

This  situation  could  occur  in  the  Marysville  project  (Figure  2.2) 
if  inflow  and  outflow  were  the  same  set  of  selective  withdrawal  intakes. 
The  effect  of  the  stratification  is  to  suppress  vertical  motions  so  that 
the  jet  must  eventually  collapse  and  entrainment  cease. 

Studies  of  the  problem  have  been  reported  by  Zu luaga- Angel ,  et  al . 
(1972)  and  Darden,  et  al.  (1975).  A  sketch  of  a  typical  How  distribu¬ 
tion  found  by  them  is  shown  in  Figure  4.9.  An  internal  hydraulic  jump 
forms  shortly  after  the  jet  enters  the  tank;  most  of  the  entrainment 
takes  place  in  the  hydraulic  jump.  The  entering  jet  pushes  ahead  of 
it  a  layer  of  reservoir  water  distinct  from  the  jet  itself.  As  a  result, 
forward  and  reverse  currents  and  rotors  were  set  up  as  shown.  Down¬ 
stream  of  the  hydraulic  jump,  the  turbulence  collapses  and  the  flow 
became  laminar  and  stable. 

The  results  of  these  experiments  enable  entrainment  to  be  estimated. 
Zuluaga-Angel ,  et  al.  (1972)  expressed  their  results  in  terms  of  the  total 
amount  of  fluid,  V  ,  entrained  in  a  time  t  after  the  discharge  was 
started.  Their  results  were  expressed  as: 

*  3/4 

E  =  Kt(t*>  1  (4.42) 

*  • k 

where  E  is  a  normalized  entrainment  parameter,  t  is  a  scaled  time, 
and  Kj  Is  a  constant. 
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(A. 43) 


,  ,2/3  11/18 

ve  v 
3/2  2/9 


qo  g 

and 

t*  -  t(eg)1/2 

(4.44) 

where 

l_^a 

p0  dy 

(4.45) 

and  v  is  the  kinematic  viscosity.  Substitution  of  Eqs.  4.44  and  4.45 
into  4.42  yields 


V 

e 


3/4,  ,-7/24  3/2  2/9  -11/18 

t  (eg)  q  g  v 


(4.46) 


To  obtain  the  rate  of  entrained  flow,  q^,  Eq.  4.46  is  differentiated  with 
respect  to  time,  yielding 


1.  -  Sr  ■  I  ‘-1M  <^)"7/24  ,03/2  62/V11/18  (A. 47, 


Finally,  a  combination  of  Eqs.  4.47  and  3.2  yields 


„  3  -1/4  ,  ,-7/24  1/2  2/9  ,-11/18 

E  -  j  Kx  t  (eg)  qQ  g  v 


(4.48) 


Zuluaga-Angel ,  et  al.  (1972)  found  two  values  for  ,  depending  on  the 
value  of  a  source  Frounde  number,  Ff ,  where 


F1 


For  210  <_  Ff  <_  985 ,  =  560 

and  for  61  <_F*  ^105,  =  250 


(4.49) 


(4.50) 


The  different  values  of  depend  on  the  strength  of  the  internal 
hydraulic  jump.  The  higher  value  of  corresponds  to  a  strong  jump. 
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the  lower  value  to  a  weak  jump. 

The  difficulty  of  estimating  the  source  of  the  entrained  fluid 
is  illustrated  by  the  complexity  of  the  circulation  patterns  induced 
(Figure  4.9).  More  flow  patterns  are  shown  in  Darden,  et  al.  (1975), 
from  which  crude  estimates  of  the  sources  of  the  entrained  flow  can 
be  made. 

Caution  should  be  exercised  in  applying  these  results  to  a  proto¬ 
type  situation.  A  dimensional  analysis  of  the  problem  yields  more 
parameters  for  the  entrainment  than  given  in  Zuluaga-Angel ,  et  al.  (1972). 
Also,  the  small  scale  and  laminar  flow  in  their  experiments,  contrasted 
with  the  turbulent  flows  in  the  prototype,  render  their  results 
questionable . 

No  studies  were  found  dealing  with  inflows  into  a  reservoir  of 
arbitrary  stratification. 

4 . 4  Buoyant  Boundary  Jet 

This  situation  is  case  C,  Figure  4.1.  The  buoyancy  forces  are 
directed  towards  the  boundary;  the  flow  is  therefore  either  a  warm  surface 
jet  or  a  cold  bottom  jet.  As  the  buoyancy  forces  oppose  the  vertical 
turbulent  fluctuations,  the  jet  collapses  at  some  point  downstream,  so 
that  E  is  always  finite.  As  this  flow  situation  occurs  frequently  in 
the  discharge  of  heated  water  from  power  plants,  many  studies  were  found. 
These  studies,  however,  concerned  only  the  discharge  into  unstratified 
ambients;  no  studies  of  discharges  into  stratified  ambients  were  found. 
4.4.1  Round  discharge,  unstratified  ambient 

The  studies  in  this  category  were  made  of  discharges  from  rec¬ 


tangular  channels.  If  the  discharges  of  interest  are  from 


similar  rectangular  channels,  the  results  can  often  be  used  directly, 
without  recourse  to  an  equivalent  semicircular  discharge  structure. 
Baddour  and  Chu  (1978)  present  a  listing  of  the  more  important 
studies  of  this  phenomenon. 

The  general  features  of  this  flow  are  shown  in  Figure  4.10.  The 

discharge  velocity  is  u,  and  the  channel  depth  and  width  h  and  2b  , 

o  o 

respectively.  Near  the  source,  the  jet  entrains  ambient  fluid  with 
the  width  and  depth  of  the  jet  growing  similar  to  a  nonbuoyant  jet. 

The  effect  of  the  jet  buoyancy,  however,  is  to  suppress  the  vertical 
turbulent  fluctuations,  so  that  the  vertical  growth  rate  becomes  less 
than  that  of  a  nonbuoyant  jet.  Eventually,  the  turbulent  flow  ceases 
entirely,  and  the  jet  collapses  into  a  thinner  layer  which  spreads 
due  to  its  buoyancy.  In  the  buoyant  spreading  layer,  the  amount  of 
entrainment  is  negligible.  Thus,  the  total  flow  entrained  by  the  jet 
is  finite  and  occurs  within  a  finite  distance  from  the  source. 

The  length  of  the  turbulent  entraining  region  can  be  estimated  from 
the  experimental  results  of  Baddour  and  Chu  (1978).  They  find  the  layer 
to  collapse  in  a  longitudinal  distance  x^,  where 


x  *  5L 
c  s 


(4.51) 


is  a  length  scale  defined  by 


,3/4 


L 


s  B  1/2 
o 


(4.52) 


where  M  is  the  kinematic  flow  force  (momentum  flux  plus  excess  hydro¬ 
static  pressure  force): 


uQ  +  Z  gfA 
o  o&  o 


4y 


(4.53) 


Bq  the  source  buoyancy  flux,  and  Qq  is  the  source  volume  flux: 

Q  =  uA  (4. 54) 

o  o 

and  Bq  *  grQo  (4.55) 

where  Aq  is  the  channel  area,  Zq  is  the  centroid  of  the  discharge,  and 
g'  is  defined  by  Eq.  4.8. 

The  amount  of  entrainment  can  be  estimated  by  the  theoretical 

analysis  of  Stolzenbach  and  Harleman  (1973),  who  considered  discharge 

from  a  rectangular  channel.  They  express  entrainment  in  terms  of  a 

maximum  dilution  D  : 

m 

Qc 

Dm  =  7T  (4.56) 

m  % 

where  D  =  f(F  ,  h  /b  )  (4.57) 

m  r’  o  o 


is  a  Froude  number  defined  by 


F 

r 


(4.58) 


Stolzenbach  and  Harleman  find  solutions  for  Eq.  4.57  by  integration  of 
the  continuity,  momentum,  and  energy  equations  and  assumption  of 
similarity  profiles.  Their  results  are  shown  in  Figure  4.11.  Combination 
of  Eqs.  3.1,  3.2,  and  4.56  yields 


E  =  D  -  1  (4.59) 

m 

Thus  E  can  be  estimated  from  Fig.  4,11  and  Eq.  4.59.  It  is  cautioned, 
however,  that  the  predictions  of  Figure  4.11  are  unconfirmed  by  experiment. 
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Figure  4.11.  Calculated  maximum  dilution,  Dm,  for  buoyant 
surface  jet  of  finite  width  (from  Stolzenbach  and  Harleman 

1973) 
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Furthermore,  the  velocity  distributions  were  found  by  Wiuff  (1978)  not 
to  be  similar,  as  assumed  by  Stolzenbach  and  Harleman.  The  reliability 
of  estimates  based  on  Figure  A. 11  are  therefore  questionable.  Stefan 
and  Shanmugham  (1975)  also  give  some  values  of  entrainment  for  x/hQ  =  10. 
These  results  cannot  be  used  here,  however,  as  entrainment  is  not 
complete  according  to  Baddour  and  Chu's  criterion. 

The  density  of  the  entrained  flow  is  indifferent  to  the  location 
of  the  entrainment  in  this  uniform  density  case. 

4.4.2  Round  discharge,  stratified  ambient 

No  experimental  or  theoretical  studies  of  this  situation  were  found. 

It  would  be  expected,  however,  that  the  amount  of  entrainment  would  not 
be  significantly  different  from  the  unstratified  case.  As  the  effect  of 
stratification  is  to  suppress  vertical  motions,  most  of  the  entrainment 
would  come  from  the  levels  of  the  discharge.  Entrainment  can  thus  be 
estimated  as  for  the  unstratified  ambient  (Section  4.4.1). 

4.4.3  Slot  discharge,  unstratified  ambient 

Many  studies  of  this  problem  have  been  made;  some  of  them  are  listed 
in  Chu  and  Vanvari  (1976). 

The  basic  properties  of  the  flow  were  analyzed  by  Koh  (1971)  and  are 
are  shown  in  Figure  4.12.  Koh  found  that  several  different  flow  regimes 
were  possible,  depending  on  the  relative  importance  of  source  buoyancy 
and  momentum  fluxes,  interfacial  shear,  and  surface  heat  exchange. 

Following  the  initial  jet-mixing  region,  an  internal  hydraulic  jump  can 
occur,  downstream  of  which  entrainment  is  negligible.  Even  if  the  jump 
does  not  occur,  the  effect  of  the  buoyancy  forces  is  to  eventually  suppress 
vertical  turbulence  and  entrainment .  Thus,  the  amount  of  entrainment  is 
always  finite. 


53 


1  rtf 

1  ENTRAINMENT 


INTERFACIAL 

SHEAR 


HYDRAULIC 

JUMP 


Figure  4.12.  General  features  of  two-dimensional  buoyant 

surface  jet 
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An  upper  limit  to  the  amount  of  entrainment  occurring  in  pumped- 
storage  flows  can  be  found  by  assuming  that  no  internal  jump  occurs. 
The  flow  is  then  of  the  maximum  entraining  type,  as  discussed  by 
Wilkinson  (1971).  Wilkinson  expresses  the  amount  of  entrainment  in 
this  case  as: 


where  F  is  defined  by  Eq. 
yields 


2F  +  1 


3F 


4/3 


4.22. 


Comparison  of  Eqs. 


(4.60) 


3.1,  3.2,  and  4.60 


E 


2F  +1 


3F 


4/3 


1 


(4.61) 


The  density  of  the  entrained  flow  is  indifferent  to  its  source  in  this 
uniform  density  case. 

4.4.4  Slot  discharge,  stratified  ambient 

No  experimental  or  theoretical  studies  of  this  situation  were 
found.  The  effect  of  stratification  would  be  expected  to  constrain  the 
entrainment  to  come  from  the  layers  immediately  below  the  discharge. 

The  amount  of  entrainment  would  probably  not  be  reduced  considerably 
from  the  unstratified  case. 

4. 5  Nonbuoyant  Boundary  Jet 

This  situation  is  case  d,  Figure  4.1.  As  there  are  no  buoyancy 
forces  exerted  on  the  jet,  the  free  surface  or  the  bottom  can  be  regarded 
as  planes  of  symmetry,  and  the  solution  techniques  are  the  same  as  those 
presented  in  Section  4.3,  with  the  amount  of  entrainment  reduced  by  one  half. 
4,5.1  Round  discharge,  unstratified  ambient 
See  Section  4,3.1. 
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4.5.2  Round  discharge,  stratified  ambient 
See  Section  4.3.2. 

4.5.3  Slot  discharge ?  unstratified  ambient 
See  Section  4.3.3. 

4.5.4  Slot  discharge,  stratified  ambient 
See  Section  4.4,4. 


5.  SUMMARY  AND  CONCLUSIONS 


The  objective  of  the  present  study  was  to  identify  mathematical 
descriptions  of  jet  entrainment  which  offer  promise  in  improving  math¬ 
ematical  modeling  of  pumped-storage  reservoir  inflows.  It  was  found  that 
a  wide  variety  of  inflow  types  could  exist,  whose  entrainment  is  influenced 
to  varying  degrees  by  source  geometry,  boundary  proximity,  source  buoyancy 
and  momentum  fluxes,  and  ambient  stratification.  Because  of  the  inter¬ 
action  of  these  many  variables,  no  one  model  exists  which  will  predict 
entrainment  in  all  possible  situations. 

The  flows  are  classified  according  to  a  primary  and  secondary  system 
discussed  in  Section  4.1.  The  primary  classification  groups  the  flows 
according  to  source  buoyancy  and  boundary  proximity.  The  four  flow  types 
are  shown  in  Figure  4.1.  The  secondary  classification  divides  each  of 
these  flows  into  those  discharging  into  uniform  or  stratified  ambient 
water  and  those  discharging  from  a  slot  or  round  orifice.  Thus,  a  total 
of  16  possible  flow  types  exist;  each  of  these  flows  is  discussed  individually 
in  Sections  4.2  through  4.5. 

Some  of  the  flow  types  have  been  thoroughly  studied,  and  some  have 
not  been  studied  at  all.  For  each  flow  type,  the  important  features  are 
discussed  and  a  candidate  method  for  predicting  entrainment  presented.  In 
general,  not  all  of  the  studies  are  cited,  nor  are  detailed  discussions 
presented.  Potential  methods  of  prediction  are  presented  either  as  analytical 
solutions,  graphical  solutions,  or  mathematical  models.  The  analytical 
and  graphical  solutions  are  given  in  the  appropriate  section;  where  can¬ 
didate  mathematical  models  are  available,  references  to  their  sourc«  ,irr 
cited. 
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The  reader  is  cautioned  that  often  only  fairly  crude  estimates  of 
entrainment  can  be  made.  This  is  partly  due  to  the  errors  involved  in 
approximating  complex  source  geometries  by  slot  or  round  orifices.  Also 
not  all  of  the  idealized  flows  used  for  entrainment  prediction  have  been 
thoroughly  studied.  As  discussed  in  the  text,  the  reliability  of  the 
predictions  for  some  of  these  flows  is  questionable.  Some  of  the  flows 
discussed  may  have  wide  applicability  to  the  prediction  of  pumped^ 
storage  effects;  these  flows  should  be  investigated  in  greater  detail, 
and  better  mathematical  models  of  them  developed. 
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APPENDIX  A:  NOTATION 


A  Channel  area 

o 

Buoyancy  flux  per  unit  length;  also  half  channel  width 
B  Width  of  slot  orifice 

B^  Source  buoyancy  flux  of  round  discharge 
Experimental  constants 
D  Diameter  of  round  orifice 

D  Maximum  dilution  of  three-dimensional  surface  buoyant  jet 

m 

E  Entrainment  coefficient 

E*  Normalized  entrainment  parameter 

F,  Fr ,  F!  Froude  number 

g  Acceleration  due  to  gravity 

gf  Modified  gravitational  acceleration 

h  Thickness  of  established  density  current  (Figure  4.2) 

h  Depth  of  channel 

o  r 

Experimental  constant 
L  Length  scale 

m^  Momentum  flux  per  unit  length;  also  dimensionless  momentum  flu 

M  Kinematic  flow  force 

M  Source  momentum  flux  of  round  discharge 

o 

Volume  flux  per  unit  length 

0,0  0  Volume  flow  rates  (see  Figure  3.2) 

<c  e,  o 

S  Dilution 

S  ,  S  Centerline  dilution 

m  o 

S  Dilution  at  terminal  rise  height 

t  Time 


A1 


t* 


Scaled  time 


T  Stratification  parameter 

u  Velocity  of  pumpback  or  generation  flow  at  source 

V  Volume  of  pumpback  or  generation  flow 

,  '  Volume  of  entrained  water 

r 

y,  z  Coordinates  (Figure  3.2) 

Length  of  entraining  flow  region 

v  Terminal  rise  height 

Z  Centroid  of  discharge 

o 

e  Stratification  parameter 

p  Dimensionless  volume  flux  parameter 

p  Terminal  volume  flux  parameter 

v  Kinematic  viscosity 

£  Dimensionless  terminal  rise  height 

p  Ambient  water  densitv 

a 

p  Ambient  water  density  of  level  of  pumpback  or  generation  flow 

and  initial  density  of  efflux  in  Figure  3.3 

p  Density  of  pumpback  or  generation  flow 
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